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Abstract. The production of the lightest Higgs in the decays of neutralinos is studied using scans of the 
phenomenological MSSM. We find that scenarios where the decay x° — > hxi is the dominant channel are 
severely constrained. However, the X2,3 ~~ ► hXi decays remain important and are highly complementary to 
multi-lepton final states in the LHC searches. The perspectives for the LHC analyses at 8 and 14 TeV as 
well as the reach of an e + e~ collider with y/s = 0.5, 1, 3 TeV are discussed. 
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1 Introduction 



The discovery of a Higgs-like boson by the ATLAS and 
CMS experiments [3H] at the CERN LHC opens up an era 
of precision studies of its production and decay properties. 
In particular, establishing if the newly discovered particle 
is the Standard Model (SM) Higgs boson or the manifesta- 
tion of an extended Higgs sector is a key question. If the 
observed Higgs-like particle is the lightest Higgs boson, 
h, of a Supersymmetric extension of the SM (SUSY), it 
could be significantly produced also in the decay chains of 
supersymmetric particles, in particular those of neutrali- 
nos, Xj ~~ ► Xih- There have been several phenomenolog- 
ical studies of h production in neutralino decays in vari- 
ous constrained supersymmetric models [5141 lj. including 
some detailed assessments of the LHC potential pj |F2"] . 
This paper discusses the regions of the MSSM parame- 
ters where these decays are relevant and the perspectives 
for their observation and study at the LHC and a future 
lepton collider in the framework of the phenomenological 
minimal supersymmetric extension of the SM (pMSSM) 
with 19 free parameters. This provides us with sufficient 
freedom to the masses and couplings to explore the su- 
persymmetric parameter space in a largely unbiased way. 
The study of the pMSSM parameter space by high statis- 
tics, flat scans of its parameters is briefly described in sec- 
tion 2. Results from measurements and searches at LEP 
and LEP2, flavour physics, dark matter experiments pro- 
vide already powerful constraints. These are discussed in 
section 3 together with the implications of the LHC results 
in the search for SUSY particles and in the first determi- 



nation of the properties of the newly discovered boson. 
Section 4 presents the regions of the viable pMSSM pa- 
rameter space probed by the contribution of the h boson in 
neutralino decays, while section 5 discusses the expected 
event rates, experimental signatures and relevant kinemat- 
ics of these processes at the LHC and a lepton collider, for 
various centre of mass energies. Section 6 has the conclu- 
sions. 



2 MSSM Scans and Tools 

The study of two-body neutralino decays into a h boson 
is based on the analysis of a large sample of MSSM points 
obtained through a flat scan of the pMSSM parameters, 
which are varied in an uncorrelated way within the follow- 
ing ranges: 



1 < tan/3 < 60, 
50 GeV <M A <2 TeV , 
-10 TeV < A f < 10 TeV, 



50 GeV < Ms 

— JL 



Mf , M 3 < 3.5 TeV . 



50 GeV < Mi,M 2 , \fi\ < 3.0 TeV 



(1) 



to generate a total of 1 x 10 8 pMSSM points. We also 
perform a scan of the CMSSM, in order to contrast the 
results for the pMSSM with those of a highly constrained 
model. For this we vary the CMSSM parameters within 
the following ranges: 



1 < tan P < 60 , 
- 10 TeV < A < 10 TeV , 

< M , Mi/2 < 3.5 TeV. (2) 

In all scans, we select the generated points which ful- 
fil a set of constraints derived from flavour physics and 
lower energy searches at LEP2 and the Tevatron, as sum- 
marised below and discussed in more details in Ref. [T5] . 
to which we refer also for further details on the pMSSM 
scans. These selected points are referred to as "accepted 
points" in the following. The pMSSM scans employ a num- 
ber of programs and software tools. Only those most rel- 
evant to this study, are mentioned here, while further 
details are given in [13]. SUSY mass spectra are gen- 
erated SOFTSUSY 3.2.3 [14] . Decay branching fractions 
are calculated using SDecay [TSIITB]. Higgs decay branch- 
ing fractions are calculated with HDECAY (5.0) [T7]. The 
flavour observables and dark matter relic density are cal- 
culated with Superlso Relic v3.2 [15] . Cross sections 
are computed with Pythia 6 . 424 [TDJ , also used for event 
generation, while the NLO k-factors are evaluated using 
Prospino v2 [2D]. The physics object response for the 
LHC analyses is simulated using the Delphes fast sim- 
ulation package [ST], tuned on the ATLAS detector per- 
formance and validated for the ATLAS multi-lepton neu- 
tralino and chargino analyses, as discussed below. 

3 Constraints 

We apply constraints from flavour physics, dark matter 
and SUSY searches at LEP, Tevatron and LHC. These 
have been discussed in details in Ref. [T3]. In particular, 
we consider the decay B s — > which can receive 

significant SUSY contributions at large values of tan/3. 
Here, we adopt the latest result from LHCb with the mea- 
surement of a branching fraction of BR(£>° — > n + = 
(3.2±i;|) x 10~ 9 with 2.2 ftr 1 of data [22], from which 
we derive the constraints at 90% C.L. after accounting for 
theoretical uncertainties. Dark matter also provides signif- 
icant bounds for this analysis. We impose the constraint 
on the neutralino relic density of 10~ 4 < Q x h 2 < 0.155, 
derived from the WMAP-7 result [23J , accounting for the- 
oretical and cosmological uncertainties and allowing the 
neutralino be responsible for only part of the observed 
dark matter. In addition, we impose the new limit on dark 
matter direct detection for spin-independent x-nucleon 
scattering recently obtained from 225 live days of the 
XENON- 100 data [23]. 

Searches for SUSY particles at the LHC in channels 
with missing Et [25, 2DJ have already provided a number 
of constraints relevant to this study, by excluding part of 
the pMSSM phase space, corresponding to gluino masses 
below ~ 600 GeV to 1 TeV and scalar quarks of the first 
two generations with masses below ~ 400 to 800 GeV. 
These constraints are included here using the same analy- 
sis discussed in Ref. [T3J, extended to an integrated lumi- 
nosity of 12 fb _1 . In addition, we also include the results of 



the searches for events with large missing transverse mo- 
mentum and two 6-tagged jets, conducted by ATLAS [27j . 
Finally, we apply the constraints in the Higgs sector by 
selecting only points having a mass of the lightest Higgs 
boson in the range 123 < Mh < 129 GeV and compatible 
with the constraints on the heavier MSSM Higgs bosons 
in the channel H/A -> r+r" [251 [2D] , 

4 Neutralino decays to h: MSSM parameter 
regions 

In the MSSM the couplings of the Higgs bosons to charginos 
and neutralinos are maximal for higgsino-gaugino mixed 
states, while the gauge boson couplings to neutralinos are 
maximal for higgsino-like states. The h couplings to neu- 
tralino are suppressed by powers of |/i|/M2 or Ma/l/il in 
the gaugino-like and higgsino-like regions, respectively. 
The branching fraction for a set of our accepted pMSSM 
points is shown as a function of \fi\/M\ and |/i|/M 2 in 
Figure [JJ 




Fig. 1. The branching fraction Xa ~> hxi for accepted pMSSM 
points as a function of the ratio \/i\/Mi (left) and |/i|/Ma 
(right). 

In the case where the charginos and neutralinos are 
gaugino-like (i.e. when the higgsino mass parameter is 
much larger than the gaugino mass parameter, \/j,\ M2) 
or higgsino-like (i.e. in the opposite situation <C M2), 
this results into the dominance of the decays of the heav- 
ier charginos and neutralinos into the lighter states and 
Higgs bosons, over the same decays with gauge boson fi- 
nal states OE]- This is also the case for the so-called "little 
cascade" in the gaugino region M2 > Mi ^> \fi\ where the 
branching fraction for the decay of the X2 the LSP 
neutralino and the lighter h boson x® ~> ^Xi, is m general 
larger than that for the decay X2 ~~ ^ -^Xij wnen kinemati- 
cally accessible in the two-body channel^] Thus, the decay 
X2 ~~ * hxl is important in these regions, unless sleptons 
are light and the additional channels where x® decays into 

1 In the three-body process, the h has to be virtual and the 
rate is suppressed by the small couplings of the h state to light 
fermions. 
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a lepton and its super-partner are open. These decays are 
most important if the neutralino is gaugino-like, since the 
coupling to the higgsino component is suppressed due to 
the small lepton mass. 

It is important to observe that for small values of 
— M2I, where the contribution of the decay into h 
is enhanced, the mass difference M v o — M^o is small, typ- 
ically of the order of 50 GeV or less. In this regime, the 
production cross sections for X2X1 an d X3X1 i n PP colli- 
sions are comparable. Due to the nature of the \2 an d X3 
states, the branching fractions into h and Z of these two 
neutralinos are complementary, i.e. the yield into h for the 
X2 is approximately equal to that into Z of the X3 states, 
and viceversa, as shown in Figure [2j This highlights the 




Fig. 2. The branching fraction X3 — > hxi as a function of that 
for X2 hxi (left) and xi — ► ^X? as a function of that for 
X2 — > Zxl (right) for accepted pMSSM points, showing the 
complementarity of the decays of xi an d X2 to h and Z. 

complementarity between the search for X2 an d X3 with 
decays into h and Z and the interest in pursuing both of 
them in the LHC analyses. 

We study the regions explorable by the X23 ~> h-Xi 
decay using our set of accepted pMSSM points. These de- 
pend only on a subset of the pMSSM parameters, mostly 
Mi, M2, j«, M2 — M(eL,Ft)- In order to characterise the 
occurrence of the various decays in different regions of 
the parameter space, we first study the fraction of the 
accepted pMSSM points into each bin in the [pi, Mi] and 
[//, M 2 ] planes where hxi, Z\1 an d the sum of££ and fr is 
the dominant two-body neutralino decay mode. We define 
the dominant decay mode as that having a branching frac- 
tion larger than both 0.2 and any of the other two-body 
decay channels. Results are shown for X2 m Figures [3] and 
|4j where we also show the region where the sum of branch- 
ing fractions for the three-body decays exceeds both 0.2 
and any of the individual two-body processes. While the 
values for the fraction of points obtained depend on the 
parameter ranges adopted for our scans given in Eq.Q, 
the picture emerging from FigureH] is of general validity. 

As we can see from Figure pT the /1X2X1 couplings 
are larger than the ZX2X1 whenever Mi > and the 
Xi is mostly bino. The region where « Mi is ex- 
cluded as the decay x§ ~~ ► hxi 1S n °t allowed kinematically. 
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Fig. 3. Fraction of accepted pMSSM points in the hi — Mi] 
MSSM plane having x% h X°i (top) Z (upper), U. or fr 
(lower) and three-body (bottom) as dominant decay. 

Very small values of |/x| are excluded from the searches 
for charginos at LEP2, which imply that \n\ > 90 GeV. 
The X2 ~~ * %x1 mode is dominant for low values of /i 
where the LSP is higgsino-like and, thus, there is a strong 
coupling of the Z boson to X1X1 states. There is a line 
Mz < M x o — M x o < Mh in which the decay channel x° — ^ 

hxi 1S kinematically close and the decay X2 ~> Zxl be- 
comes dominant, provided that the sleptons are not light 
enough for the ££ channel to contribute. Finally, at low 
H values, Xi an d X2 are higgsino-like, but the splitting is 
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Fig. 4. Same as Figure^ for the [fj, - M 2 ] MSSM plane. 



typically small, so the two-body decays through h or Z 
are forbidden and the three-body modes dominate. 

Through the decays h — > bb, h — > t + t~, Z — > bb, 
Z —> t + t~ , these processes lead to the final states bbl 
+ E™ ssins and ttI + £™ issing at the LHC, where the 
X2 3 is produced in association with a chargino decaying 
xf — > W ± Xi — > ^ ± ^Xi- An alternative, intriguing sce- 
nario leading to the same ttI + _E™ Issmg final state arises 
if the fi is light. In this case, the decay x° ~~ > — * TT Xi 
may be enhanced. 



5 Neutralino decays to h: pMSSM to 
CMSSM comparison 

Wh production from X\ X% decays at the LHC has been re- 
cently discussed in the context of the CMSSM in Ref. [30], 
where the authors find that the decay x% ~> ^Xi is dom- 
inant for typical choices of the model parameters. In our 
analysis we find indeed that X2 ~~ ^ ^Xi is the dominant 
decay, with branching fraction in excess of 0.20, for 71% 
of the accepted CMSSM points before applying the relic 
density constraint. This is no longer the case when im- 
posing the full set of constraints. Of the accepted points 
with X2 ~^ hXi as the dominant decay, only 10~ 4 ful- 
fil the constraint < Q x h 2 < 0.155 and 8xl0~ 5 the 
tighter requirement of 0.068 < Sl x h 2 < 0.155 set by the 
WMAP bound. In fact, the region where the x° ~~ hXi 
is dominant corresponds to values of neutralino relic den- 
sity, f2 x h 2 , which largely exceed the upper limit of the 
WMAP result |23j . as shown in Figure [5] This is due 
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Fig. 5. Branching fraction for X2 — ► hxi as a function of the 
neutralino relic density, fl x h 2 , for pMSSM (left) and CMSSM 
(right) points. 

to the peculiar nature of the CMSSM, which has large 
/j. values, with 3> Mi making the neutralinos to be 
gaugino-like, as observed above. The general prevalence 
of the x,2 ~^ hXi decay, observed in the CMSSM before 
the fi x h 2 constraint, cannot be generalised to the MSSM. 
But in turn the greater flexibility of the pMSSM allows to 
reconcile the relic density constraint to large coupling of 
the h to the neutralinos (see Figure[5|. For comparison, in 
the pMSSM the decay into is dominant in only 10% 
of the accepted points before the relic density constraint. 
But now 28% of these fulfil also the loose relic density con- 
strain of 10~ 4 < Q x h 2 < 0.155 and still 2.2% make also 
the tighter requirement of 0.068 < Q x h 2 < 0.155. It is 
therefore important to reconsider the apparent dominance 
of the xl -> X°h decay in the CMSSM, when examining 
the SUSY phenomenology in the context of more general 
models, or even in the CMSSM itself after imposing the 
dark matter relic density constraint. Still, the process of h 
production in neutralino decays is of crucial importance in 
specific regions of the parameter space, which are largely 
complementary to those yielding decays into sleptons or Z 
bosons, probed by the multi-lepton analyses at the LHC. 
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6 Neutralino decays to h: Experimental 
Signatures 

In the following we discuss the perspectives for detecting 
and studying these decays at the LHC and an e + e~ linear 
collider (LC). We consider for the LHC an integrated lu- 
minosity of 25 fb- 1 at 8 TeV and 300 fbr 1 at 14 TeV. For 
the LC we consider 500 fb _1 of integrated luminosity at 
= 0.5 TeV and 1 atr 1 at 1, 1.5 and 3 TeV. For each of 
these scenarios, we combine the decay branching fractions 
with the relevant neutralino production cross sections to 
study the number of expected events. For the case of the 
LHC, we also characterise the kinematics of the detectable 
final states and study a possible event selection to isolate 
neutralino decays into hxx (and Z\\). We compare these 
results to multi-lepton final states. 

6.1 LHC 

As a first step, we combine the results on the branch- 
ing fractions discussed in the previous section to the ex- 
pected relevant cross sections for the associated chargino 
neutralino production. Through the decay xt — > W ± x±, 
W ± —ti^v, this process produces a highp^ lepton, which 
is very valuable for the event trigger and subsequent selec- 
tion. The cross sections for chargino neutralino production 
at 8 and 14 TeV, followed by X2.3 ~* ^Xi! h ~> bb, tt and 
xf — > W ± Xi, W ± — tr^v are given in Figures [6] and [jj 
respectively. 

Then, we study the fraction of accepted pMSSM points 
in bins of the Mi, M2 and [i parameters, which yield more 
than 10 signal events for x° — * ^X? an d X3 ^Xi with 
25 ft)" 1 at 8 TeV and 300 fb _1 at 14 TeV. Results are sum- 
marised in Figures [8] and [9j respectively. Again the values 
for the fractions of points given in these figures depend on 
the the ranges adopted in the scan at Eq.Q. However, the 
location of the regions highlighted by the large occurence 
of points with enough signal events are of general valid- 
ity. Due to the factor of ~3 increase of the pp — > X1X2 3 
production cross section and the factor of ~10 increase 
in the assumed statistics, the coverage of the parameter 
space expands significantly when moving from 25 fb _1 at 
8 TeV to 300 fb" 1 at 14 TeV. Still, the 8 TeV data provide 
interesting limits, when comparing to the other channels 
sensitive to SUSY weak production at the LHC. Results 
from searches for neutralino and chargino production have 
been reported for the two- and tri-lepton channels by both 
ATLAS and CMS [5IH55] . These channels are sensitive 
to the contribution of either light sleptons in the decays 
X° -> £t; £ -> iXi, or Z bosons in the decay x° ^Xil 
Z — > 11. In the unconstrained MSSM, slepton masses can 
be pushed well above the X2 mass, so that the only remain- 
ing dominant two-body x® decays are either hx® or Zxl- 
This motivates us to pursue the study of the challenging 
h X i channel, through the bb£ + £™ issing and t+t~£ + 
^missing t p i gj eSj wn ich complement the Zx\ channel. 
In addition, the decays X2 ~^ ^ T are 01 special importance 
in the region of light f 1 , which is highlighted by a possible 



enhancement of h 77 yields in the ATLAS and CMS 
data. This process leads to the same t + t~£ + i?™ lssmg 
final state as the decays /iXii h —> tt and Zxl', Z tt. 
If a signal is observed in the bbl + E™ 1SSU1S and t + t~£ + 

E T lns modes it would become essential to attempt to 
identify the contributing channels, possibly through a fit 
to the bb invariant mass and the tt transverse mass. 




Fig. 6. Cross section for X1X2 ( u PP er ) and xtxi (lower) 
production vs. M x o and M x o, respectively, with ^2,3 hxi', 

h -> bb, tt and xt -» ^xi; W /± -> in 8 TeV pp colli- 
sions for a set of accepted pMSSM points. 



In order to compare the coverage of the space of the 
relevant pMSSM parameters for the Higgs and the multi- 
lepton channels, we compare the results obtained above 
to the fraction of the accepted pMSSM points for which a 
signal can be obtained in the two- and three-lepton chan- 
nels across the [/i — Mi] and [fi — M 2 ] planes at 8 and 
14 TeV, based on the signal selection of the ATLAS anal- 
yses. Accepted pMSSM scan points fulfilling the selection 
criteria discussed above and having an inclusive chargino 
and neutralino production cross section larger than 0.01 fb 
at 8 TeV collision energy are studied in details. A set of 20k 
SUSY events with inclusive neutralino and chargino pro- 
duction is generated for each point at both 8 and 14 TeV. 
Simulated events are passed through a fast parametric 
simulation using the Delphes package which simulates the 
physics objects used in the subsequent analysis. The se- 
lection criteria of the two- and three-lepton + missing Et 
ATLAS analyses [3TJ [35] are applied. The observability of 
the pMSSM point is assessed based on the correspond- 
ing CI value, where the number of expected signal events 
is obtained from the result of the selection on the gener- 
ated events rescaled to a luminosity of 25 fb^ 1 and that 
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Fig. 8. Fraction of accepted pMSSM points in the [fj, — Mi] 
and \p, — M2] planes giving at least 20 events in the hxii^Xi 
channel for 25 fb" 1 of data at 8 TeV. 



of background events is taken from that estimated for the 
ATLAS analyses, rescaled again to the assumed luminos- 
ity. The results of this simulation is validated by compar- 
ing the 95 % C.L. exclusion contours obtained in the M 2 
- /1 plane with those expected for the ATLAS analyses, 
under the same assumptions. These contours agree within 




H (GeV) 

Fig. 10. Fraction of accepted pMSSM points in the \p — Mi] 
and [fi — M2] planes observable with 25 fb" 1 of data at 8 TeV 
in the three-lepton analysis. 

20-30 %. Results for 25 ftr 1 at 8 TeV and 300 fb" 1 at 
14 TeV are given in Figure [10] 

Results for searches in the bb£ + i?™ Il!Slns channels have 
not been reported yet by the LHC experiments. In order 
to obtain a first assessment of the observability of the bb£ 
+ _E™ lsslng final state over the pMSSM parameter space, 
we perform a simple analysis on signal events generated 



G 




Fig. 11. Distribution of the average reconstructed event trans- 
verse missing energy, _E™ laslng , (left) and bb di-jet energy, Ebb, 
(right) in signal events at 14 TeV from the accepted pMSSM 
points. 



at our pMSSM points with a h boson produced in a X2 3 
decay, based on the Delphes fast simulation response. In 
this analysis, jets are reconstructed using the anti-kt algo- 
rithm |Mj with a cut 0.4. Leptons are selected requiring 
Pt > 10 GeV and pseudo-rapidity below 2.40 for elec- 
trons and 2.47 for muons and to be isolated. We select 
events with only one lepton and two jets with Ex > 25 or 
45 GeV. We further require the reconstructed jet to be as- 
sociated to a generated b hadron. These selection criteria 
are inspired to those adopted in the WH Higgs analysis 
with H — > bb, leading to a similar final state, but with 
smaller values of E T llsslns . Here, we do not investigate the 
backgrounds, for which the ti and W + jets processes 
are expected to give the dominant contributions. In our 
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Fig. 12. Fraction of accepted pMSSM points in the [fi — Mi] 
giving more than 20 selected signal events at 14 TeV in the 
bb£ + _B™ ls&ms analysis described in the text for an _B™ lsslng cut 
of 50 GeV (upper) and 100 GeV (lower plot). 



signal analysis, we require significant missing transverse 
energy, i5™ Isslng in the event, due to the two escaping x?s. 
The average value of the transverse missing energy and 
of the average energy of the reconstructed bb system in 
hxl, h — » bb for signal events from the pMSSM 



A2,3 



points which pass the cuts above, is shown in Figure 11 



We perform the analysis for values of the cuts on the b 
jets p T of 25 and 45 GeV and on E™ issing of 50, 100 and 
150 GeV. The efficiencies at 14 TeV for b jet p T 25 GeV 
are typically in the range 15 t o 3 0% depending on the re- 



quirement on £?™ lssmg . Figure 
pMSSM points in the M x 
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shows the fraction of the 



parameter space where the 
analysis selects more than 20 signal events. In the case the 
X2 an d X3 states are close in mass, they would be both 
produced and they may decay into hxl and Zx\, respec- 
tively, as discussed above. In this case the bb£ + _B™ Issing 
final state will receive contribution from both h — » bb and 
Z — » bb, which can be resolved at least on a statistical ba- 
sis from an analysis of the invariant mass of the bb system. 



6.2 e + e Linear Collider 

An e + e~ collider of sufficient energy and luminosity is 
very well suited for reconstructing two-body decays of 
pair-produced neutralinos, e + e~ — > XiX% including those 
decaying into h and Z bosons. In particular, since the en- 
ergy of the produced Xi is known, apart from beam-related 
effects, the energy distribution of the reconstructed boson 



can be used to precisely determine the Xi mass. The de 



Zx\ was studied as a benchmark reactions at 



cay xl 

■\fs = 0.5 TeV using full simulation and reconstruction as 
part of the ILC Lois [351 EE]- I n these analyses the X2 
mass of 217 GeV was determined with a statistical rela- 
tive accuracy of 1%. The decays X2 ~> ^Xi an d X2 ~^ Zx\ 
have been studied for CLIC at s/s = 3.0 TeV [37]. The 
X2 —> ]}Xi decay, followed by h — > bb, leads to the signa- 
ture bbbb+ missing energy final state. The X2 decay can 
be isolated in an inclusive SUSY samples, where its mass 
and h yield are determined with a statistical accuracy of 
4% and 5%, respectively. 

The pair production cross sections multiplied by the 
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X2 3 — ► hxi branching fractions are shown in Figure 
a function of the neutralino mass, for the four y/s energy 
values (0.5, 1.0, 1.5 and 3.0 TeV) chosen for this study. 

The production of neutralino pairs in e + e~ collisions 
is mediated by s-channel Z boson exchanges as well as t 
and u-channel left- and right-handed e exchange. If the 
sleptons are assumed to be very heavy or the produced 
neutralinos are higgsino-like (and thus the coupling x°ee 
is negligibly small) only the s-channel Z boson exchange 
is relevant. In fact, the neutralino cross sections are rather 
small, even in the case of X1X2 production where the phase 
space is more favourable. The only exception is when both 
Xi and X2 are higgsino-like and have maximal couplings 
to the Z boson. Also the cross sections with identical neu- 
tralinos are extremely small. Due to Fermi statistics, the 
neutralinos are produced in p- waves and therefore sup- 
pressed at threshold and, in the case of gauginos or hig- 
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Fig. 13. Product of neutralino pair production cross section 
and neutralino decay branching fraction into hxi as a function 
of the X2 (upper) and \3 (lower) masses at ^fs — 0.5, 1, 1.5 
and 3 TeV. 



pMSSM points where the sum of the neutralino masses is 
more than 100 GeV below the assumed collision energy, 
to avoid threshold effects. With an expected integrated lu- 
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Fig. 14. Neutralino pair production cross section for e + e~ — > 
X2X2 an d X2X3 a t \f& = 1-5 TeV shown as a function of the 
X2 ~~ ► h-Xi branching fraction for pMSSM points where the sum 
of the neutralino masses does not exceed 1.4 TeV. 



gsinos, the cross section is further suppressed by the cou- 
plings Za oc Nf 3 — Nf A . The cross section for mixed X1X3J 
X2X3 and X3X4 are significant, except in the case where 
one of the neutralinos is a pure gaugino, which leads to 
suppressed Zxx couplings. It is important to observe that, 
in the region of parameters where the x\ — > h\i process is 
the dominant neutralino decay, the x§ tends to decouple 
from the Z. This suppresses the e + e~ — > X2X2 P air P ro ~ 
duction, leaving instead the e + e~ — > X2X3 process as the 
only significant neutralino production mechanism, when 
kinematically available. This correlation between the pair 
production cross section and the decay branching frac- 



tions is demonstrated in Figure 14 for a sample of accepted 



Fig. 15. Fraction of accepted pMSSM points in the [fj, — Mi] 
plane giving at least 20 events in the hxihxi channel for 1 ab _1 
of e + e~ data at -Vs = 1.5 and 3 TeV. 



minosity of 0.5 - 1.0 ab _1 a high energy linear collider can 
observe these decays over a significant part of the kine- 
matically allowed region of the parameter space. This is 
quantified here by the fraction of the pMSSM points yield- 
ing at least 20 signal events for an integrated luminosity 
of 0.5 ab- 1 for = 0.5 TeV and of 1 ab" 1 for y/s = 
1, 1.5 and 3 TeV, over the [fi, Mi] and [/i, M 2 ] parameter 
space, as shown in Figure fl5| 

Finally, we study the determination of the yield of 
h bosons in an inclusive SUSY sample for the specific 
case of a MSSM model having with mass of 340 GeV, 
X2 with mass of 643 GeV decaying predominantly into 
h\i and \i with mass of 643 GeV decaying exclusively 
into W ± x° 1 and M h = 126.1 GeV, in 3 TeV e+e" colli- 
sions. An inclusive sample of SUSY events, corresponding 
to 0.5 fb _1 of integrated luminosity, is generated using 
Pythia and events are fully simulated using Geant-4 [55] 
and reconstructed following the same analysis procedure 
discussed in |37j . Events are pre-selected requiring a visi- 
ble energy 0.08-^/s < E to t < 0.6-^/s, an energy in charged 
particles larger than 150 GeV, transverse energy larger 
than 200 GeV and jet multiplicity 2< Nj ets <5. Jet clus- 
tering is performed using the Durham jet algorithm |39j , 
with y cut — 0.0025, on the reconstructed particle flow ob- 
jects of the Pandora particle flow package [ID]. Selected 
events are clustered into four jets and the di-jet invari- 
ant mass for all the three possible pairings is computed. 
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The jet pairing minimising the difference between the di- 
jet invariant masses is selected, provided the mass differ- 
ence is below 20 GeV. The fraction of background non- 
resonant events is obtained from the di-jet mass side- 
bands 20 < E j:j < 60 GeV and 140 < E j3 < 200 GeV. 
The fraction of W + W~ , Z°Z° and h°h° events in the in- 
clusive 4-jet SUSY events is extracted by a x 2 fit to this 
di-jet mass distribution. The W and Z° mass peaks are 
parametrised as Breit-Wigner functions convoluted with 
a Gaussian term describing the experimental resolution. 
The mass and width values of the Breit-Wigner functions 
are fixed to their generated values, while the total area and 
the width of the Gaussian resolution terms are left free in 
the fit. The h° peak, which has negligible natural width, 




100 120 140 161 

Di-jet Invariant Mass (GeV) 




120 140 1G0 180 
Di-jet Invariant Mass (GeV) 



Fig. 16. Di-jet invariant mass in inclusive 4-jet + missing en- 
ergy SUSY events produced in y/s=3 TeV e + e~ collisions for 
all selected events (upper) and 6-tagged events (lower). The 
result of the fit is shown by the continuous line with the in- 
dividual W, Z and h components represented by the dotted 
lines. 



is modelled as the sum of two Gaussian curves, one repre- 
senting the correctly reconstructed signal events, centred 
at the nominal Mh value, the second describing decays 
where the mass has a lower reconstructed value due to 
semi-leptonic b decays. The central value, width and frac- 
tion of events in this second Gaussian is extracted by a fit 
to a pure sample of decays into h bosons and fixed in the 
fit, while the Gaussian width of the main peak is kept free. 
The yield of h bosons, with mass 126.1 GeV, is extracted 
by a x 2 fit to the di-jet mass distribution of events in 4- 

for all di-jets and for 



jet final states, shown in Figure 16 



those passing di-jet b-tagging based on the ZVTOP algo- 
rithm [5T] implemented in the LCFIVertex package [4"2~| . 
We measure the fraction of h bosons to be 0.269±0.013 



and that of Z bosons to be 0.037±0.016, after the non- 
resonant background subtraction, which compares well to 
the original values of 0.290 and 0.025, respectively, of the 
generated events. This specific example shows the feasibil- 
ity to accurately determine the yield of h bosons in SUSY 
particle decays in the data of an e + e~ collider operating 
at sufficient energy. 



7 Conclusions 

The observation of a Higgs-like particle with mass ~ 126 GeV 
promotes its study in the decay of supersymmetric parti- 
cles. Neutralino decays offer a prime opportunity, since the 
decay x\ 3 — \ hx\ is highly complementary to the chan- 
nels leading to multi-lepton final states, such as Z and 11. 
The typical products of cross sections and decay branch- 
ing fractions to the bbl final state is of the order 0.01 - 
1 fb in 8 TeV pp collisions and a factor of ~ three larger 
at 14 TeV for M v o < 600 GeV. This should make pos- 
sible to search for this decay in the bbl + MET over a 
significant region of the parameter space at 14 TeV, com- 
plementary to that covered by searches with multi-leptons, 
provided the ti and W + jets backgrounds can be properly 
suppressed. An e + e~ linear collider of sufficient energy, 
y/s > 1 TeV, and luminosity, yielding >0.5 ab _1 of data, 
can study h production in neutralino decays for masses 
up to the kinematical limit for pair production, obtaining 
0(1%) precision on the mass and 0(5%) on the yield of 
Higgs bosons from neutralino decays in inclusive multi-jet 
SUSY events. 
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